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ABSTRACT 


If de Sitter' 8 hydrostatic equations are developed independent 
of the external potential theory , the hydrostatic geopotential 
coefficient J. occurs explicitly on the right-hand side of those 
equations, since J ^ here has to be treated as an unknown in the 
solution , it becomes rather difficult to solve the equations 
independently , regardless of which of the dynamical parameters 
associated with the earth is taken as the initial datum. Solution 
is possible , however, with the help of a boundary condition derived 
from the external potential theory which neither assumes nor discounts 
the presence of equilibrium conditions in the earth's interior. If 
a general solution is constructed on this basis , the three particular 
solutions, usually quoted in literature, stem from it in the wake 
of the appropriate assumptions. Of course, the only meaningful 
solution — of these— is that corresponding to the polar moment of 
inertia as the initial datum. It is essential that the solution be 
constructed in this way in order to demonstrate clearly the correct 
structure of the problem of hydrostatic equilibrium. 

The anomalous gravity field of the earth referred to the hydro- 
static figure is compared with that referred to the international 
reference ellipsoid. 




Introduction 


The hydrostatic theory of the earth, in its present form, was 
originally developed by Clairaut (1743), Radau (1885) simplified 
it by making an important substitution in Clairaut 1 s differential 
equation. Originally it was hoped that a knowledge of flattening f, 
dynamical flattening H, and m (defined later) would lead to some 
useful information about the distribution of density within the earth. 

The application of the theory, however, showed that with in and H 
fixed, widely different la^s of density led to almost the same 
value of flattening ( Jeffreys , 1962, page 152). This led to an impor- 
tant result: if one could assume that this theory was applicable 

to the real earth, i.e., if the earth was in hydrostatic equilibrium, 
it would be possible to compute its flattening. This would require 
a knowledge of the parameters m and H and of the distribution of 
density, to compute the geopotential coefficient J. The parameters 
m and H were quite well-determined. The distribution of density 
was not very well-known, but this theory had indicated that the use 
of a reasonable density model was adequate for the computations. An 
obvious inference was that the value of flattening, obtained on the 
basis of this theory for a reasonable density distribution, would give 
a good approximation to the actual flattening of the earth spheroid. 

Such an inference was particularly welcome at that time because a 
precise determination of the geopotential coefficient J was not 
possible to provide the value of flattening directly. Other methods 
of determining the flattening of the earth were dependent upon gravimetric 
or arc measurement data and hence had their obvious limitations. Doubts 
as to the validity of the above assumption arose, however, when 
Tisserand (1891) and Poincare (1910) derived f, 1 « 297.3 (f, denotes 

— ft ft 

hydrostatic flattening and f real flattening) on the basis of hydro- 
static theory, in contrast to the then-accepted value of f ^ = 293.5 
which A.R, Clark had obtained from arc measurements in 1880 ( Jeffreys , 
1962, p. 152). Later, Hayf ord (1909, 1910) and Helmert (1911), using 
more extensive data, obtained f ^ - 297.0, which was quite consistent 
with the value obtained from hydrostatic theory within the limits of 
accuracy. This reaffirmed the then-prevalent opinion that the method of 
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hydrostatic theory for computing the flattening of the earth was more 
accurate than any other available at that time ( de Sitter , 1924, 1938; 
Jeffreys , 1948, 1962; Jones , 1954), This interesting conclusion was 
based on the belief that since hydrostatic equilibrium probably 
existed throughout the earth’s interior except for the crust and perhaps 
the upper mantle, the assumption of hydrostatic equilibrium for the 
whole earth was not unjustified. This belief persisted until artificial 
satellites made possible the direct determination of the geopotential 
coefficient J, from which the flattening f of the best-fitting spheroid 
could be accurately computed using the external potential theory alone, 
without recourse to either gravimetric data or geodetic arc measurement 
data or to the hydrostatic equilibrium theory. The currently accepted 
value of f computed in this way is 1/298.23. Since the real flattening 
and the hydrostatic flattening could now be computed independent of 
each other, any discrepancy between the real flattening f and the 
hydrostatic flattening f^ could now be studied, and the validity of 
the claim made in the pre-satellite times that f - f^ could now be 
proved or disproved definitely. Henrik sen (1960) used de Sitter ’s 
equations and J * 1622.4 x 10 ^ to find a value of 1/300 for hydrostatic 
flattening. Jeffreys (1963), using J = 1624.17 x 10 ^ as obtained by 
King-Hele , Cook and Rees (1963), found a value of 1/299.67 ± 0.05 for the 
hydrostatic flattening, confirming the discrepancy between f and f^ 
which Henriksen (1960) and O’Keefe (1960) had initially pointed out. 

Using the improved value of the constant of precession ( Rabe , 1950), of 
moon-earth mass ratio (1/81.303) and J = 1623.969 x 10 ^ ( Kozai , 1964), 

I ( Khan , 1967) obtained 1/299.86 ± 0.05 for the hydrostatic flattening 
by employing Jeffreys ’ (1963) approach, which has both the advantages 
of speed and simplicity. In his method, Jeffreys (1963) has carried the 
analytical development to first order and then applied the corrections 
for the second-order terms by the numerical integration of the appro- 
priate expressions. Although Jeffreys ’ paper outlines the structure of 
the hydrostatic problem clearly, he does not use de Sitter ’s equations. 
There has been considerable confusion regarding the computation of 
hydrostatic flattening from de Sitter ’s hydrostatic equations which 
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carry the development of the theory to the second order analytically. The 
general belief has been that the hydrostatic flattening is and should be 
computed from the hydrostatic part of the equations alone, to the 
complete exclusion of the external potential theory. In order to 
examine the validity of this belief I re-evaluated (Khan, 1968a) the 
hydrostatic theory and made some modifications in de Sitter 's (1924) 
equations in order to reinstate the geopotential coefficient J which 
de Sitter had eliminated from the right-hand side of his equations by 
means of a relation derived from the external potential theory. This 
step is essential in order to demonstrate the correct structure of the 
problem. When de Sitter worked out his equations, the geopotential 
coefficient J was not known accurately and he could presume that there 
was no significant difference between the value of J for the real earth 
and the value of hydrostatic J. This assumption has been proven as 
incorrect ever since. However, as will be discussed later, there is 
theoretical justification for the elimination of J' — but the reason that 
we ought to avoid it at this stage is that it is so well-concealed in 
the development of de Sitter ’s equations that it tends to obscure the 
correct structure of solution of the problem of hydrostatic equilibrium. 

In this report, I first give the derivation of the modified 
equations followed by a description of the general solution. 1 then 
present a review of previous methods, both pre-satellite and post- 
satellite. 

I wish to emphasize here that the discussion given in this report 
pertains primarily to de Sitter 's development of the hydrostatic theory. 
The elimination of the geopotential coefficient J from the right-hand 
side of de Sitter *s equations is certainly valid if the correct reason 
for it is outlined clearly. However, the prevailing situation is that 
it has given rise to the belief that the hydrostatic flattening can be 
computed merely from the hydrostatic part of de Sitter 's development of 
the theory without any use whatsoever of the external potential theory 
even if the only known datum used is the polar moment of inertia of the 
real earth. The principal purpose of this report is to show that if 
de Sitter ’s development is used for computing hydrostatic flattening 
merely from a knowledge of the earth's polar or mean moment of inertia 
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(making no additional use of the satellite-determined J except in the 
computation of the polar moment of inertia), some control from the 
external potential theory must be introduced. 

Theory of Hydrostatic Equilibrium 

Development of Modified Hydrostatic Equations: 

The customary manner of representing the outer potential of a 
reference figure, symmetrical with reference to both an equatorial 
plane and the polar axis (to the order of accuracy required here) is 


GM 
V = — 

1 - J 2 

'V! 

L. 

P 2 (sin <}>) - 

H 

a 

— P 4 (sin *) 

r 


lr 1 


r 1 


'1 2 2 2 , y , <v 

+ - oj r coa <(> (1) 


where 


G « gravitational constant 
M « mass of the earth 
a e = equatorial radius of the earth 
a) = rate of rotation of the earth 

(sin <}>) and (sin <j>) * Legendre’s polynomials 
and - zonal harmonic coefficients in the spherical 
harmonic representation of earth’s potential 

The radius vector of an ellipsoid of revolution can be represented 
(to the order of accuracy required here) as 

“ a f a + a 0 P 0 + a P ) 
e K o 22 4 4' 


r 


( 2 ) 
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where 


a * 
o 










a 


4 


12 .2 

35 f 


(3) 


la the foregoing expression f is the flattening of the surface 
described by Eq. (2), 

Substituting r from Eqs. (2) and (3) into Eq. (1) and retaining 
2 

quantities to the 0(f ) only, we get the potential V on the surface of 
the ellipsoid, i.e.. 


V = V q (GM, a e , m, f, J 2 ) + (GM, a 0 , m, f, J 2 > 

6 ^ (GM, a^, m, f, J 2 , J^) 

If the ellipsoidal surface is an equipotential surface, V must be 
constant on it, and the coefficients B 2 and 6 ^ must both be zero, which, 
among other relations, gives us the following equation 

2 12 2 1 

J *= “T" f — "T" f + — - mf - -r- m (4) 

2 3 3 21 3 w 

In the derivation of Eq. (4), no assumption is made regarding the 
conditions existing in the earth’s interior. Hence the equation is 
valid whether or not hydrostatic equilibrium exists in the earth’s 
interior. Eq. (4) gives 


* 3 1 , 15 T 9 t 2 3 2 

f - j J 2 + j m + — J 2 m + g J 2 + 5 ^ m 


(5) 
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Equations (4) and (5) should relate the observed J with the real 
flattening in a non-hydrostatic case and the hydrostatic J with the 
hydrostatic flattening for a hydrostatic case. Hence Eq. (4) or (5) 
will be of the form 


G (J, f) ■« 0 


(5a) 


for the non-hydrostatic case and 


G (J h ,f h ) - 0 


(5b) 


for the hydrostatic case. 

Note that and f^ denote the hydrostatic J and the hydrostatic 
flattening, respectively. 

The real flattening of the ellipsoidal surface which best fits 
the geoidis given by Eq. (5) in the form of Eq. (5a), if J for the 
real earth is known accurately as, for example, from the regression of 
node of an artificial earth satellite. 

The condition of hydrostatic equilibrium can be represented 
analytically by the differential equation 


dp = 6dV 


where V is the sum of gravitational and rotational potential, P the pres- 
sure at any point within the body and 6 its density. This implies 
that p must be a function of V, and 6 either a constant or a function 
of V also. Since the equipotential surfaces are given by V - constant, 
p and 6 will also assume a constant value, as obvious from Eq. (3). 

Thus, the condition of hydrostatic equilibrium for the earth’s interior 
is established by stipulating that the equipotential surfaces and the 
surfaces of equal density be coincident. If 6 represents the mean 
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radius of any such surface, expressed in terms of the outer surface as 
a unit so that 8 ranges from 0 at the center to 1 at the outer surface, 
then the equation of this surface becomes 


r - e [l - j (f + | f h 2 ) P 2 (sin f) + ~ f h 2 P 4 (sin $)] (6) 


whe re 


f f 


J5 

42 



Since we are developing a second order theory and thus are not 

2 2 2 

interested in terms > 0(f ), we will treat f* = f^ and mf * = mf^, 
etc* 

The potential V at any point (r,$) within the earth is given by 
(de Sitter , 1924): 


GW 


D(1 + j n^o 3 cos 2 <fr) - | (So 2 + T o 3 ) P £ (sin ♦) 


12 -4 32 5 

+ Qo' ) P 4 (sin <f>) 


(7) 


where 


2 3 

a) a m 

e 

= - 

GM D 

W = volume of the surface of which 8 is the mean radius 
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D 


S 


T 


P 


Q 


6 

6 e 2 dB 


B 5 (f + I f h 2 ) 


/ S “ «’ + 21 £ h^ dB 
J a QP 





dS 


dB 


In these formulas 6 is the density, expressed in terms of mean 
density as a unit, and D is the mean density within the surface B, 
expressed in the same unit. For the outer surface 


D - 1 
s 


S s - 3 J < 1 + ! f h> 


P 

s 


35 

12 


J 


4 


T - 0 
s 

Q » 0 
s 


( 8 ) 


The condition that the surface, described by Eq. (6), be an 
equipotential surface, gives as one of the relations the following: 

D (f 1 + y f h 2 - i m 1 ) - j (S + T) = 2i f (m - 3T) 


(9) 
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Differentiating Eq. (9) and Introducing a new variable 
by 


d log V £ df' 
0 ” d log B = f' dB 


we get 

- 3D(1 + f f h ) - ff- 

With the help of Eq. (8), the above equation gives us a 
condition for the outer surface which can be written as 


D(1 + 


f h> ~ 


_4 

21 


m + 


n s f ' (1 + 7 £ h - 21 "> ■ 3f,<1 + 7 f h> - 5J h (1 + 3 


where n denotes the surface value of n* 
s 

Again, by differentiating Eq. (9a), we can get 
6 dB + ^ + 5n ~ + n) - 21 « * 0 


in which 


C 


1 dD 
D dB 



and 


C “ 7 | (1 + n) - 3f h (l + n ) 2 - 4f 


» defined 
(9a) 

boundary 
f h ) (10) 

(ID 

(11a) 


(Hb) 
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Radau (1885) transformed Eq. (11) to the form 


d_ 

dB 


DB 5 Vl + n = 5D B 4 F(n) 


( 12 ) 


where 


F(n) 



1 2 
10 n 



Vl + n' 


(12 a) 


F(n) has the remarkable property that it lies very near unity, 
it f s maximum departure being 0.0007* If 1 + A denotes the average 
value of this function over the range of integration, we can write 


f D B 4 dB - ~ — *• 1 1 B 5 

J Q 1 + * 


which, for the outer surface, assumes the form 




1 + A 

s 


(13 a) 


Note that n and 
s 

outer surface. 

Now in our units, 


A are the values of parameters n and A for the 
s 

we have 



4 71 
3 


a + f y 


The polar moment of inertia C can be written as 


_8 

15 


X 

* / 6 dB t® 5 (1 + 3 Vl dB = 3 * / 


+ f (c - A) 


6 B 4 dB 
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where A is the moment of inertia about the equatorial diameter. With 

C - A 

H = 

C 



and hence 


J 


h 


q H , 


we obtain another boundary condition for the outer surface in the form 


q = 3 


1 - — f +A.£ 2 
3 h 9 h 


j 

1 1 


6 B 4 dB + 4 J, 
3 h 


(14) 


As stated before, if D is the mean density within the surface 
8 expressed in terms of the mean density as a unit, we have 


3 B 

- I 


6 8 dB 


which on differentiation with respect to g gives 


ri dD = 3 1 _ & 

D dg D 


or 


6 = D |i + Ii® 

0 | 3 D dg 
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Substituting this value of 6 in Eq. (14) we get 


q = 3 


3 l h 9 


£ h 2 ) / D ( 1 + Tif) 64dB 




which, with the help of Eq. (13a) can be finally transformed into 


, 2 2 4 „ 2 2 

1 - f f. + 4 J, + f f, - f 
3 h 3 h 9 h 5 


i £ +i( M 

3 h 9 h I 1 + X 

S 


(15) 


Equation (10) can be simplified to the form 


V - 


3f ' - £ 

J 7 


m f, 


5 + 


10 

21 


, , 20 
f h + 21 m 


(16) 


From Eq. (15) one obtains 


I 33 

A 

,2 

1 

S i 



i I,. 2 T . 4.2 

1-q ~3 f h + 3 J h + 9 f h 


2 2 4 2 

— 1 - — f + — f 

5 3 h 9 h 


- 1 


(17) 


At this point, if we can assign a reasonable value to A (for 

s 

example, X = 0 as explained later), we can substitute n from 
s s 

Eq. (17) into Eq. (16) and obtain an equation in and f^ from which 
a series of values of f^ can be obtained corresponding to arbitrary 
values of J^. However, it is convenient to get an expression which 
gives f^ explicitly in terms of other parameters. This can be done 
by writing Eq. (17) as 


25 _2 ,2 

— F 9 



( 18 ) 


where 
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q 1 * l - q 


II 




. 2 I 

h I 


A i If 2 2 

2*3 n ~ 3 r h 


F - 1 + X 


( 19 ) 


It is instructive to note that A^ and A^ are both of the order 
of flattening. 

Simplifying Eq. (18), one obtains 


25 _2 ,2 

n - —r F q' 
s 4 n 


1 + 2(A 2 - A x ) 


+ | A x 2 + 3A 2 2 - 4A 


a)] - 1 ■ »o + °1 + "J 


(20) 


where 


2fF 2 ,' 2 -l 


25 -.2 ,2 

T F « ( *2 


V 


25 _2 , 

— F q' 


, * . 2 ,2 

J h • q f h + 3 ’ £ h 


and 


( 21 ) 


2| r 2 ,.2 




34. 


- 41 1 4 2 


Note that the quantity is of the order of f^, whereas r\^ 

is of the order of f, 

n 
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Using this value of n , Eq. (16) can be written as 

s 

A f h 2 + (n 0 " 3 + 6 1> f h + 5J h + 6 2 " 0 < 22 > 

where we have put 


, _ 17 _5 25 v 2 , 

A 14 42 n 0 " 3 F q q 


{ 1-T r2 ’' J h " t ” + 21 J h 


(23) 


and 


6 


2 



/ 


Note that 5 is approximately of the order of f, , whereas <5. is 
* o n 2 

of the order of f, . 

h 

Equation (22) gives the required expression for f, which, correct 

h 

to the second order of small quantities, is 


1 



- (5J h + V + 


25 A J. 2 
n 

(n 0 - 3) 2 


(5j h + 6 2> 6 l 
(n 0 - 3) 


(24) 


It is interesting to see that in this development, the expression 
for ^corresponding to the first order theory is 


f 


h 



3 - 


n 


0 


(25) 



-17- 


General Solution of the Problem of Hydrostatic Equilibrium 

To solve the problem of hydrostatic equilibrium we have to examine 
hydrostatic Eq. (24) which can be written as 


F(m, H, J h , f h ) = 0 


(26) 


In terms of the basic parameters, the above equation can be written as 


F(C, A, f^» to, a, M) *» 0 


(26a) 


If the hydrostatic equilibrium exists, the figure of the earth predicted 
from the external potential theory should be coincident with the 
hydrostatic figure (See Footnote, p. 37 ) • This gives an important boundary 
condition of the problem, i.e., if the hydrostatic equilibrium exists in 
the earth’s interior, Eqs. (5) and (24) must match at the outer boundary 
of the earth. Hence, 


F(m, H, J h , f h ) = G(J h , f h ) 


or more specifically. 


3 - n f 


(5J + 6.) 6 25 A J ‘ 

<5J. + 6.) + — ! 2 ^ i + h_ 

n Z- -» / n \ 4 


3 - n 


0 


(3 - n Q > 


( 27 ) 


■ I J h 2 + P 1 J h + F 2 


1 + 


_5 

14 


m 


where 
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Fig. 1. Graphical solution of Eq. (9) for (cj, M, a, H) * constant and (a), M, a, C) = constant. 


OIOG 


19 



Fig. 2. Graphical solution of Eq. (9) for (w, M, a, J) * constant. 


- 20 - 


v 1 . 3 2 

F 2 = 2 m + 56 m 

Other quantities on the right hand side of Eq. (24) are defined 
previously. 

Equation (27) gives the general solution of the problem of 
hydrostatic equilibrium. The particular solutions are obtained by 
examining the left-hand side of Eq. (27) in the form of Eq. (26) or 
(26a) and by properly defining the basic parameters occurring in that 
equation. If the rate of rotation u) and the mass M are chosen to be the 
same for the hydrostatic earth and the real earth, there are three 
possible solutions which correspond to the following boundary conditions 


(to. 

M, 

a, 

C) = constant 


(ll), 

M, 

a, 

H) = constant 

(28) 

(u), 

M, 

a. 

J) - constant 



Graphical solutions of the general hydrostatic Eq. (27) are given in 

Figure 1 for the first two cases and in Figure 2 for the third case. 

The plots of the function G(J^, f^) are shown by solid curves, whereas 

those of F(m, H, J^, f^) subject to the set of boundary conditions 

mentioned in Eq. (28) are shown by broken-line curves. It is obvious 

from the shape of these plots that each set has an unique intersection 

point in the region in which we are interested. Note that in Figure 1 

2 

the abcissa is scaled in terms of where J n , = *r J, . In all the 

Zn Zn J ix 

solutions, A g has been taken equal to zero. It is shown later (see 

Fig. 3) that reasonable variations in values of A do not affect the 

s 

solution critically. The rate of rotation a) is treated as constant 
via the parameter m in all the solutions given in Figures 1 and 2. 

For the case when (to* M, a, C) = constant, the hydrostatic 

flattening is f ^ - 299.75 ± 0.05 and the hydrostatic J is J « 1607.49 

-6 ^ “ 
x 10 . This value is about the same as reported by Henriksen (1960) , 
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x$ 

Fig. 3. Influence of A on the value of f, . 

s h 


O’Keefe (1960), and in a previous paper by me ( Khan , 1967). However, 
as much as one can gather from the post-satellite literature, the method 
reported previously ( O’Keefe in footnote to Munk & MacDonald , 1960) 
certainly does not profess ( Khan , 1968a, b) to take care of the boundary 
condition noted above. The fact that the two methods give identical 
numerical results can be traced to de Sitter ’s specific elimination of 
the geopotential coefficient J from the hydrostatic equations in order 
to obtain his equations for the hydrostatic theory. This elimination 
remained undetected ( Khan , 1968a, b) when initial results on hydrostatic 


flattening were reported and led to the erroneous belief that the then- 



- 22 - 


reported value of hydrostatic flattening was obtained from a solution 
of the hydrostatic equation alone. The net effect of this elimination 

is, in fact, the same as allowing for the above-noted boundary condition. 

However, it should be clear now that if the above value of f, is pre- 

h 

ferred, the method reported in this report is the correct way to obtain 

it. As usually believed, for this model the hydrostatic flattening is 
smaller than the real flattening. 

The merit of this solution lies in the fact that the polar moment 
of inertia of the hydrostatic earth is equal to that of the real earth, 
as obtained from observational data on the geopotential coefficient J 
and the dynamical flattening H. Consequently, it is possible to avoid 
certain dynamical complications (mentioned later) which arise when the 
moment of inertia of the hydrostatic model is taken different from that 
of the real earth. 

For the case when (a), M, a, H) = constant, the hydrostatic flattening 
is f^ - 297.29 ± 0.05 and the corresponding hydrostatic J is = 1623.225 

x 10 This value of hydrostatic flattening is very near the flattening 

obtained in pre-satellite times. The pre-satellite method essentially 
consisted of predicting J from the hydrostatic theory and using it in 
Eq. (5) to compute hydrostatic flattening which was then assumed to give 
the best approximation to the real flattening. Thus, the pre-satellite 
method will tend to give similar results as the solution proposed here, 
if the only data used in the solution is the dynamical flattening H. 

Note that the precise determination of the geopotential coefficient J has 
not changed the method of solution. It has merely made possible the de- 
termination of the earth's moment of inertia which provides a better initial 
datum than the dynamical flattening H. 

However, the polar moment of inertia of this hydrostatic model is 
greater than the real earth — and this introduces some dynamical 
complications. It would imply a change in the radial stratification of 
the earth, such as would result because of the equatorial bulge of the 
real earth being more compressed and consequently the real earth having 
a higher density gradient than the hydrostatic state would require. 

This creates the problem of suggesting some reasonable physical 
phenomenon responsible for such a process. Some increase in the polar 
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moment of inertia could possibly be accounted for by the fact that 
when the earth readjusts to the equilibrium shape defined by the above 
model, there will be an increase in the polar moment of inertia because 
of the expansion of the equatorial bulge to conform to the new figure. 
Approximate calculations show, however, that this factor can account 
only for a small fraction of the total variation required by this model. 
This is the hydrostatic model of pre-satellite times when the geo- 
potential coefficient J for the real earth was not precisely known and 
hence, the moment of inertia of the real earth could not be determined. 

For the hydrostatic model (m, M, a, J) - constant, the solution of 
the equation G(J^, f^) = 0 will obviously give a constant value of 
flattening, as can be seen from Figure 2. The hydrostatic flattening 

for this model is f 1 = 298.29 ± 0.05 and the corresponding hydrostatic 

n 

value of H is H = 3260.50 x 10 . However, the polar moment of inertia 

of the hydrostatic model is greater than that of the real earth. The 
dynamical problems encountered in this case are again of a similar 
nature to those enumerated for the second hydrostatic model. 

The results of the above three solutions are summarized in Table 1. 


, TABLE 1. Possible Hydrostatic Figures Obtained 
from the General Solution _a, and Data for the Real Earth b 





a_ 






Hydrostatic 

Model 

J , x 10 6 
n 

H x 10 6 

C/Ma 2 

A/Ma 2 

f h 

-1 


v 1 - '- 1 

(u), M, a, H) 

1635.225 

3273. 64 

0.33300851 

0.33191836 

297.29 

+ 

0.05* 

-0.96 

(w, M, a, C) 

1607.49 

3240.43 

0.33071598 

0.32964432 

299 . 75 

± 

0.05* 

+1.50 

(oj, M, a, J) 

1623.969 

3260.50 

0.33204876 

0.33096611 

298.29 

± 

0.05* 

0 




b 







J x 10 6 

H x 10 6 

C/Ma 2 

A/Ma 2 

f" 

-1 




1623 . 969 t 

3273.64** 

0.33071598 

0.32963333 

298.25 

t 

0.05 



* Based on m » 0.00344980 ( Khan . 1967; Jeffreys . 1964) 
t Kozai (1964) 


** *than (1967) 
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Choiee of A 
s 

In the above discussion the value of A is taken as zero. As 

s 

pointed out earlier, the function F(n) always lies very near unity, 

-4 

its maximum deviation being of the order of 10 . A g which is the value 

of A for the outer surface, is given by Bullard (1948) as 

X g = (1.6 ± 1.8)10~ 4 

This estimate, however, is based on the density distribution 

suggested by Bullen (1940, 1942). Jeffreys (1963), using a simplified 

-4 

density model, finds A = 1.3 x 10 and points out that if A « 0 

S y S 

instead, the resulting f^ is greater by 6 x 10 only. Khan (1968a) 

showed that slight variations in the value of A do not affect the 

s 

value of to any great extent, as is apparent from Figure 3. 

It seems legitimate, therefore, to take A g = 0 for initial calculations. 

Note that Henriks en (1960) also took A g - 0. 

Table 2 summarizes some of the more important values of A and n 

s s 

as obtained by different investigators. 

TABLE 2. Comparison of Some Hydrostatic Theory Parameters 




X 

s 

n s 

de Sitter 

(1924) 

(4.4 ± 1.5) ID -4 
-4 

0.5589 

Bullard 

(1948) 

(1.6 ± 1.8)10 
-4 

0.565 

Jeffreys 

(1963) 

1.3 x 10 4 

0.5587 

Khan 

(1967) 

— 

0.5869 
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Minimum Strength of the Earth 

The stress differences arising because of the departure of the 
earth from hydrostatic equilibrium are given by Jeffreys (1943, 

1963) . On the supposition that the stresses are supported by strength 
(1) down to the core, or (2) down to a depth of 0.1 of the earth ? s 
radius, the strength S needed to support the inequality is given as 
follows: 

For the hydrostatic model ( w, M, a, C) = constant: 

Case 1 

S - 4.3 x AJ^ x 10^ dynes/ cm^ = 4.7 x 10^ dynes/cm^ 

Case 2 

S = 7.9 x AJ^ x 10 ^ dynes/ cm^ - 8.7 x 10^ dynes/ cnt 


Gravity Field Referred to the Equilibrium Figure 

The first hydrostatic model (obtained from the earth's polar 
moment of inertia as the datum) is used as reference for computing the 
anomalous gravity field of the earth from Kozai 9 s (1964) zonal and 
Gaposhkin 8 s (1966) tessera! harmonic coefficients. If the anomalous 
gravity field is to be used for any studies regarding the earth 9 s 
crust and mantle, it must be computed with reference to the equili- 
brium figure ( O'Keefe and Kaula , 1963; 0 s Keefe ^ 1965; Fischer , 1967) 
of the earth, because such a figure is a figure of zero stress, and 
departures from it will, inter alia, be indicative of the hydrostatic 
stresses existing in the earth f s crust and mantle. Figure 4 shows 
the gravity anomalies referred to an ellipsoid with flattening 1/299.75. 
The gravity anomalies referred to the international reference ellipsoid 
are shown in Figure 5. It is obvious from a comparison of Figure 4 
with Figure 5 how the picture of the anomalous gravity field is a 
function of the model adopted as the reference « This is also evident 
from some of the discussions given by O’Keefe and Kaula (1963) and 
Fischer (1967). In any case it is clear that the satellite-determined 
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Fig. 5. Gravity anomalies in milligals referred to the international reference ellipsoid. Computed 
using Kozai's (1964) zonal and Gaposhkin's (1966) tesseral harmonic coefficients to P OQ . 
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gravity anomalies referred to any of the reference models do not 


exceed about 43 milligals for an 8th ham 

ionic representatioi 

n of the 

gravity field* as seen from Figures 4 anc 

1 5. The most pron 

ounce d 

gravity anomaly in Figure 5 is negative i 

md occurs in the Ii 

ndian 

Ocean just to the south of Ceylon. The n 

lagnitude of this a; 

nomaly 

reduces by about 11 milligals to about 32 

\ milligals (Fig. 4 

) when the 


equilibrium figure is used as the reference* However, in that case 
the positive gravity anomaly over the islands of New Guinea and Borneo 
becomes accentuated by an almost equal amount and becomes the most 
pronounced feature. Also the negative gravity anomaly located to the 
east of Zapadno and Sibirskaya in USSR becomes more pronounced by 
about 11 milligals. The well-pronounced positive anomaly to the south- 
southwest of Iceland in Figure 5 is less pronounced in Figure 4. The 
two negative anomalies flanking the southern tip of North America are 
equally observable in both the representations, while the positive 
anomaly over and around Peru is much more pronounced in Figure 4. The 
negative anomaly over the Hudson Bay area is decidedly more pronounced 
in Figure 4. Several other contrasting features of interest can be 
pointed out from a study of the two gravity field representations. It 
is interesting to compare the variance of the anomalous gravity field 
with respect to the various reference figures: 

Variance of the anomalous gravity field with respect to the 

2 

equilibrium figure - 145 mgal . 

Variance of the anomalous gravity field with respect to the inter- 

2 

national reference ellipsoid = 136 mgal . 

Figures 6 and 7 show the geoidal undulations with respect to 
the international reference ellipsoid and the equilibrium figure, 
respectively. Figure 8 is the astrogeoid computed with reference to a 
flattening of 1/299.67 and is taken from Fischer (1967). A comparison 
of Figure 7 with Figure 8 shows some very interesting features. The 
geoid depression at latitude 60° north of Tomsk in USSR is coincident 
in location in Figures 7 and 8 but it does not seem to match very 
closely in magnitude. However, there is no astrogeodetic data farther 
north of latitude 60° in the location of this particular depression and 
it is difficult to say whether or not this depression will deepen farther 




Fig. 7. Geoidal undulations (in meters) with respect to the equilibrium figure (f ~ * 299*75) * Computed 

using Kozai's (1964) zonal and Gaposhkin f s (1966) tesseral harmonic coefficients to P 
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Fig. 8. Astrogeoid (in meters) referred to the equilibrium figure (f, = 299.67) . (From Fischer, 

1967.) h 
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north . The geoidal low over Hudson Bay is equally well represented in 
both figures, A comparison of the geoidal high over the northwest and 
west coast of South America in Figure 7 with its counterpart in Figure 
8 is strikingly noteworthy. Over Africa the two geoidal representations 
do not agree too well, but it should be noted that astrogeodetic data 

for the continent of Africa are very sparse. The two well-pronounced 
features of the satellite geoid, i.e., the geoidal low in the Indian 
Ocean and the geoidal high over the Solomon Islands area, cannot be 
compared because practically no astrogeodetic data are available for 
those areas. The agreement in the major features of the two geoidal 
representations computed from two independent sources is, among other 
things, indicative of the appropriateness of the figure of reference 
to which the satellite and the astrogeodetic geoid are referred. 


Review of the Previous Methods 

Below we give a very brief review of the previous methods, examining 
them in perspective in the light of the general solution. 

Pre -Satellite method : 

If we eliminate from the right-hand side of Eq. (15) and (16) 
with the help of Eq. (4), we obtain 


, 1 , 1 ^ 2 _ 2 - 

4 3 9 h 21 h 


2,. 2 .4,2. 

5 (1 3 f h + 9 f h 5 


V 1 + 


(29) 


1 + A 


a 5 , '10 2 4 - 2 6 f 

V = 2”- 2f + 2l m + 7 f h “ 7 “ f h 


(30) 


The second order terms in Eq. (29), i.e., the terms containing 

2 3 

and m f^, ultimately become 0(f ) when they are multiplied by f on 

substitution in Eq. (30) and drop out. Their inclusion is desirable, 

however, to show exactly where the different terms become negligible. 

With the exception of the second-order terms, Eq. (29) and (30) are the 

same as de Sitter’ s equations (21) and (22) in his 1924 paper. 



-33- 


The method adopted by de Sitter is discussed in detail in numerous 

papers ( de Sitter , 1924, 1938; Bullard , 1948; Khan , 1968a). H is taken 

as the initial datum and an attempt is made to find a value of 

which would be compatible with the selected value of H (both J, and H being 

■ 2 

functions of m) . This is done by estimating a quantity q = 3/2 C/Ma^ 

(J^ = qH) from a knowledge of the internal density distribution of the 
earth ( de Sitter , 1924; Bullard , 1948). Once is determined, Eq. (5) 
will give the corresponding hydrostatic flattening. A serious disadvan- 
tage of the solution obtained by this method is that the moment of inertia 
of the hydrostatic model, constructed in this way, is greater than that 
of the real earth. This creates some dynamical complications as discussed 
before. However, it must be appreciated that d£ Sit ter *s whole effort 
was really directed to devise a method which would give the best approxi- 
mation to the real flattening of the earth, not necessarily the hydrostatic 
flattening, in the context of the information available in his time. 

Other post-satellite methods : 

Jeffreys (1963) has given an excellent numerical method based on a 
simplified density model. He computes the various hydrostatic parameters 
from the first-order theory and evaluates the second-order correction 
terms by the numerical evaluation of the appropriate integrals. The 
method reported in this paper could be really regarded as a counterpart 
of Jeffreys 1 (1963) method with the exception that I have employed 
de Sit ter Vs development of the hydrostatic theory. 

In the previous applications of de Sitter ! s development to compute 
the hydrostatic flattening, however, it is usually claimed that using 
satellite-determined J and dynamical flattening H, and hence 
knowing the polar moment of inertia of the earth, f^ should be computed 
from the hydrostatic equations alone without using any controls from 
the external potential theory. If this is accepted, the use of de 
Sitter T s (1924) equations is automatically ruled out because these 
equations are derived with the help of external potential theory. 
Consequently, the modified equations ( Khan , 1968a) should be used, but 


in that event, the equations cannot be solved because of the explicit 
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TABLE 3. Comparison of Hydrostatic Flattening Values* 


Pre -Satellite Method 

(Based on the dynamical flattening H as the initial datum) 



de Sitter 

( 1924) 

296.92 

± 0.136 

de Sitter 

(1938) 

296.753 

± 0.086 

Bullard 

(1948) 

297.338 

± 0.050 

Jeffreys 

(1952) 

297.299 

± 0.071 

Jeffreys 

(1963)t 

296.75 

± 0.05 


Results Obtained After the Satellite ^determination of J 
f^ - hydrostatic flattening; f - 1/298*25 ± 0*05 


(a) 

Based on 

the polar moment of inertia as 

the initial datum 




V 1 

'h- 1 - '■ 

Henrik, sen 

(1960) 

300.0 


+1.75 

0 'Keefe 

(1960) 

299.8 


+1.55 

Jeffreys 

(1963) 

299.67 

± 0.05 

+1.42 

Khan 

(1967) 

299.86 

± 0.05 

+1.61 


(b) From a solution of the modified equations alone, 
using polar moment of inertia and the satellite-determined J 

296.70 * 0.05 ** -1.55 

297.04 t Q.OStt -1.21 

*See also Table 1. 

tUsing the pre-satellite approach. 

**Baaed on m - 0.00344980 ( Khan , 1967a; Jeffreys , 1963), H - 0,00327364 ( Khan , 
1967), and J 2 « 0.001082645 (Kozai, 1964). 
ttBased on m * 0.00344992 ( Henrikaen , I960), H * 0,00327070, and 
J 2 - 0.00108270. 
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appearance of on the right-hand side of the equations and because 
of the necessity of treating this quantity as an unknown in the solution. 
If these modified hydrostatic equations are solved with the help of 
satellite-determined J (i.e., J - J^) one will obtain the results 
shown in the last section of Table 2. The first part of Table 2 gives 
the results obtained from different investigators using de Sitter ’s 
equations in a way which is pertinent to the discussion given in 
this report. 

Table 3 lists the results obtained from the modified hydrostatic 
equations ( Khan , 1968a) along with the results of other investigators 
who obtained their values using de Sitter s equations. 

The important work of Ledersteger (1967) is not discussed because 
the discussion given here does not pertain to his method. 

Summary and Conclusions 

Since the hydrostatic geopotential coefficient J, appearing 
explicitly on the right-hand side of the modified hydrostatic 
equations, must be treated as an unknown in the solution of the 
hydrostatic equilibrium problem, these equations become difficult 
to solve all by themselves and one has to look for an additional 
boundary condition. This boundary condition is inherent in the definition 
of the hydrostatic equilibrium and is derived from the external 
potential theory which neither assumes nor discounts the existence 
of hydrostatic equilibrium in the earth f s interior. It requires that 
the equilibrium figure of the earth coincide exactly with that pre- 
dicted from the external potential theory, if hydrostatic equilibrium 
exists in the earth’s interior, and is stated in terms of Eq. (5)* 

The solution obtained with the help of this boundary condition turns 
out to be sufficiently general so that the three most frequently 
mentioned particular solutions in literature can be obtained from this 
by merely defining the appropriate initial datum. The only geophysically 
meaningful model, of course, remains the one in which the polar or 
mean moment of inertia is held equal to that of the real earth, 
calculated from the satellite-determined J and the dynamical flattening 
H computed via the constant of precession of the real earth. For 
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this hydrostatic model the flattening is = 299*75. The solution 

is significant in that it demonstrates the correct structure of the 

problem of hydrostatic equilibrium of the earth. Taking the equilibrium 

figure as reference, the minimum strength of the earth required to 

support the hydrostatic stresses (considering the inequality due to 

7 2 

second zonal harmonic only) is 4.7 x 10 dynes/ cm in case the stresses 

7 2 

are supported by strength down to the core, and 8.7 x 10 dynes/cm 
in case the stresses are supported by strength down to a depth of 0.1 
of the earth’s radius. 

The anomalous gravity field of the earth with respect to the 
equilibrium figure and the international reference ellipsoid is shown 
in Figures 4 and 5. The equilibrium figure provides the best reference 
for computing the anomalous field because such a field would also 

reflect the hydrostatic stresses which become very important in 
geophysical studies on a regional scale. A comparison of Figures 4 
and 5 shows the extent to which the anomalous gravity field depends on 
the reference figure used. Figures 6 and 7 show geoidal undulations 
with respect to the international reference ellipsoid and the equi- 
librium figure. Figure 8 shows the astrogeoid referred to the equi- 
librium figure. The agreement between the various extremums in Figures 
7 and 8 is noteworthy and reflects the appropriateness of the reference 
figure to which both geoids are referred. 
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(Footnote for p. 17) 

It was brought to the author's attention that this was also recognized 
by Caput o (1965), who compared the hydrostatic flattening obtained from 
several different assumptions. However, Caput o did not point out the 
problem mentioned in this report, i.e., the incompleteness of the 
system of equations (15 and (16), if is treated as an unknown, and, 
hence, does not deal with the problem considered in this report. 



BIBLIOGRAPHY 


Bullard , E , C . , The figure of the earth, M. N . R .A .S . , Geophys . Sup p 1 . 
5, 186-192, 1948. 

Sullen, K. E . , The problem of the earth’s density variation, Bull . 

Seism. Soc. Am ., 30, 235-250, 1940, 

Bullen, K. E . , The density distribution of the earth’s central core. 
Bull. Seism. Soc. Am ., 32, 19-29 , 1942. 

Caputo, M., The minimum strength of the earth, J. Geophys . Res ., 70, 
955-963, 1965. 

Clair aut , A. C., Theorie de la Figure de la Terre , Paris, 1743, 
de Sitter, W., On the flattening and the constitution of the earth. 
Bull. Astr. Insts . Netherlands , 2, 55, 97-108, 1924. 
de Sitter, W., On the system of astronomical constants, Bull. Astr . 

Insts . Netherlands , 8, 307, 213-229, 1938. 

Fischer, I., Deviations of the geoid from an equilibrium figure, 

Army Map Service, Washington, D. C., 1967. 

Gaposhkin, E. M., A dynamical solution for the tesseral harmonics of 
the geopotential and for station coordinates (abstract) , Trans . 
Am. Geophys. Union , 47, 47, 1966. 

Hayford, J. F., The figure of the earth and isostasy, U ,S.C.&,G.S . , 

1909. 

Hayford, J. F., Supplementary investigations in 1909, U..S.& S.G.S ., 

1910. 

Helmert, F. R . , S . B. P r eu so , Akad . Wi ss . , 10-19 , 1911. 

Henriksen, S. W. , The hydrostatic flattening of the earth, Ann. Int . 
Geophys. Year , 12, 197-198, 1960. 

Jeffreys, H., The stress differences in the earth’s shell, M,»N .R.A.S 
Geophys. Suppl. 5, 71-89, 1943. 

Jeffreys, H. , The figure of the earth and the moon, M.N ,R.A.S . , 
Geophys. Suppl. 5, 219-247, 1948. 

Jeffreys, H., The Earth , Cambridge University Press, Cambridge, 1952 
1962. 

Jeffreys, H. , On the hydrostatic theory of the figure of the earth, 
Geophys. J ., 8, 196-202, 1963. 


38 - 



”39- 


Jones, H. S., Dimensions and rotation, Chapt. I, In G. P. Kuiper (ed.) 

The E ar th as a Planet , The University of Chicago Press, Chicago, 

1954. 

Khan, M. A., Some parameters of a hydrostatic earth, Trans. Am. Geophys . 
Union, 48, 1, 56, 1967. 

Khan, M, A., A re-evaluation of the theory for the hydrostatic 

figure of the earth, J . Geophys . Res . , 73, 5335-5343, 1968a. 

Khan, M. A., General solution of the problem of hydrostatic equilibrium 
of the earth, (in press), 1968b. 

King-Hele, D. G. , G. E. Cook, and J. M. Rees, Determination of the even 
harmonics in the earth’s gravitational potential, Geophys. J ., 8, 

119, 1963. 

Kozai, Y., New determination of zonal harmonic coefficients of the earth’s 
gravitational potential. Smith. As trophy s. Qbsy. Sp. Rept . 

No. 165, 1964. 

Ledersteger, K. , Critical notes on mean functions derived from orbital 
perturbations of artificial satellites. In G. Veis (ed.) , The Use 
of Artificial Satellites for Geodesy, vol. II , National Technical 
University, Athens, Greece, 1967. 

Munk, W. H. and G. J. F. MacDonald, Con tinen tali ty and the gravitational 
field of the earth, J. Geophys. Res ., 65, 7, 2169-2172, 1960. 

O’Keefe, J. A., Determination of the earth’s gravitational field, In_ 

H. Kallmann (ed.). Space Research I . 448-457, 1960. 

O’Keefe, J. A., Letter to the Symposium on the determination of the figure 
of the earth, Studia Geophys ica et Geodaetica , 2, 9, 1965. 

O’Keefe, J. A. and W. M. Kaula, Stress differences and the reference 
ellipsoid. Science , 142, 3590, 1963. 

Poincare, H., Bull. Astronomique France , 27, 321-356, 1910. 

Rabe, E . , The p recessional constant from Eros, Astronom. J ., 55, 181, 1950. 

Radau, R., Compte Rendus De Science, C, R. Acad. Sci. Paris , 100, 972-974, 
1885. 

Tisserand, F . , Mecanique Celeste , 4 Vo Is . , Paris, 1889-1896. 



